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ABSTRACT 

We re-examine critically the arguments raised against the theory that Ultra 
High Energy Cosmic Rays observed at Earth are produced in Gamma Ray Bursts. 
These include the limitations to the highest energy attainable by protons around 
the bursts' shocks, the spectral slope at the highest energies, the total energy 
released in non-thermal particles, the occurrence of doublets and triplets in the 
data reported by AGASA. We show that, to within the uncertainties in our 
current knowledge of GRBs, none of these objections is really fatal to the scenario. 
In particular, we show that the total energy budget of GRBs easily accounts for 
the energy injection rate necessary to account for UHECRs as observed at Earth. 
We also compute the expected particle spectrum at Earth, showing that it fits 
the HiRes and AGASA data to within statistical uncertainties. We consider the 
existence of multiplets in AGASA' data. To this end, we present a Langevin-like 
treatment for the motion of a charged particle in the IGM magnetic field, which 
allows us to estimate both the average and the rms timedelay for particles of given 
energy; we discuss when particles of identical energies reach the Earth in bunches, 
or spread over the rms timedelay, showing that multiplets pose no problem for 
an explosive model for the sources of UHECRs. We compare our model with 
a scenario where the particles are accelerated at internal shocks, underlining 
differences and advantages of particle acceleration at external shocks. 

Subject headings: shock waves - cosmic rays - gamma ray: bursts 



1. Introduction 

We suggested some time ago (Vietri 1995) that Gamma Ray Bursts provided the ideal 
accelerating sites for UHECRs, supporting this view with a computation of the highest en- 
ergy attainable by protons accelerated around the fireball's external shock (which was found 
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to exceed lO'^^eV), and of the average energy deposited by bursts per unit time and volume 
in gamma ray photons (which came close to that inferred for UHECRs, if these have a cosmo- 
logical origin); assuming comparable efficiencies in the production of gamma ray photons and 
non-thermal particles, we argued that this provided an excellent explanation, without fitting 
any parameter, of the UHECRs' fiux levels at Earth. Since then, our knowledge of bursts 
has increased dramatically, and the fireball theory (Rees and Meszaros 1992) has received 
a spectacular confirmation with the accurate description of the properties of afterglows (see 
Piran 1999, 2000 for reviews). Independently of this, the recently reported data from HiRes 
(Abu-Zayyad et al, 2002) appear once again to support the view that the origin of the 
Ultra High Energy Cosmic Rays (UHECRs), at least as observed at Earth, is to be found 
among astrophysical sources. These two facts together make it worthwhile to re-examine 
the objections that were raised against the above-mentioned scenario, in the wake of its first 
appearance. Also, an alternative model for the production of UHECRs in GRBs has been 
proposed (Waxman 1995), and thus it seems pertinent to discuss some differences between 
these two models, which will be subject to different observational tests. 

The plan of the paper is as follows: in Section 2, we address the criticisms levied against 
the model in the past years which concern the acceleration mechanism. In Section 3, we deal 
with global energetics. We discuss in Section 4 two caveats on the previous discussion on 
energetics. Doublets and triplet in the data reported by AGASA, are discussed in Section 
5; we compare our work to that of other authors in Section 6, then we summarize our 
conclusions. 



2. Acceleration 

The details of the model for the acceleration of non-thermal particles remain essentially 
those outlined in Vietri (1995), and in Dermer (2002a): particles are accelerated around 
the fireball's external shock, which is highly relativistic for most of the relevant time. The 
major reason for this is that the basic fireball model employed (Meszaros, Laguna, and Rees, 
1993, where a detailed model for the external shock evolution is presented, together with the 
physical characteristics of the plasma behind the shock) has proved itself able to explain the 
characteristics of the afterglow (Piran 1999, 2000), and it thus appears in need of no major 
revision. The model was used mostly in determining the largest energy which could bounce 
off the ejecta shell, without crossing it unscathed (and also, in showing that no radiative 
losses where incurred into, by the highest energy particles). We derived for the highest 
energy 
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which differs from the original (Eq. 34 in Vietri 1995) because a more modern initial Lorentz 

factor r = 100 has been used, the true energy release E and beaming angle 6 have been 
introduced. The above equation defines the largest energy which can be turned back in a 
typical wind which produces GRBs, as seen in the upstream frame (the laboratory); the 
very important question of why we should see particles leaving the shock region through 
the upstream frame rather than the more customary downstream frame will be addressed 
shortly; let us just remark that particles leaving through the downstream region will be seen 
in the laboratory frame as having an energy roughly smaller than the above. 

The average length for deflecting back a non-thermal particle has been taken as ^ 40 
times its Larnior radius, but the computation holds only for F ^ 1. It was shown in Vietri 
(1995) that losses due to all radiative processes are negligible. Furthermore, no particle 
loss may occur through the sides of the ejecta (which would be conceivable, because the 
emission is beamed along a cone of semiopening angle 9); in fact, in the comoving frame, 
the ejecta have a (radial) width given by r/4F, where r is the shock's distance from the 
explosion site, while the ejecta width in the transverse direction is ~ r6, and, since 1/F < ^ 
(Panaitescu and Kumar 2002) it will be easier for the particles to leave through the bottom 
end than through the sides. Recently, Panaitescu and Kumar (2002) have obtained values 
for the circumburst's density n, for the departures of the magnetic field from its equipartition 
value 63, for the total energy E release and beam opening angle 9, for a variety of bursts's 
afterglows for which sufficient time and frequency coverage were available; using the data 
in their Tables 2 and 3, we compile our Table 1, illustrating the largest energies achievable 
within this scenario, which are conveniently large with respect to observational data. 
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TABLE 1. Largest energy achievable by protons accelerated at external shocks, and 

particle spectral index. 



GRB 


£^ma./102° eV 


970508 


7.2 


980519 


38.0 


990123 


5740. 


990510 


161. 


991208 


1.80 


991216 


18. 


000301c 


0.5 


000418 


0.2 


000926 


1.4 


010222 


1.3 
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A critique against this argument was levied by Gallant and Achterberg (1999), who 
pointed out that the limiting factor for the highest energy achievable is due to the ability of 
the insterstellar medium to deflect high-energy particles toward the shock, and not the vice 
versa, which is what was computed above. Their analysis showed that this highest energy 
was 

E^ax ~ W . (2) 

This argument however assumes a definite value for the magnetic field around the burst, 
which is very model-dependent: there is no compelling reason to assume for B the value of 
the average ISM. In particular, two of the leading models for GRBs, the binary neutron star 
merger (Narayan, Paczynski and Piran 1992), and the SupraNova (Vietri and Stella 1998) 
are clearly at odds with this value. In fact, these two models make very similar predictions, 
because in both cases the GRB goes off inside a pulsar wind bubble (PWB, Konigl and 
Granot 2002; Inoue, Guetta and Pacini 2001, Guetta and Granot 2002, but the same, bare 
proposal was made also in Vietri 1995, and in Gallant and Achterberg 1999), filled with an 
electron/positron gas, and its associated magnetic field. In either case, a pulsar-type dipolar 
magnetic field {B = 10^^ — 10^^ G) has a light cylinder at ric ~ 10^ cm, and decreases like 
B oc r^^ (at least in its toroidal component), from there on, as dictated by standard pulsar 
electrodynamics (Contopoulos, Kazanas and Fcndt, 1999). This results in fields of order 
i? = 0.1 — 10 G at the radii appropriate for afterglows (r = 10^^ — 10^'^ cm). These estimates 
are in agreement with determinations of the magnetic field of the (much slower rotating!) 
Vela bubble (Helfand, Gotthelf and Halpern 2001, but see Arons 2002 for a different estimate 
of the magnetic field) . 

In order to adapt Eq. 2 to our needs, we repeat Gallant and Achterberg's computations 
for B oc r~^. We consider a particle leaving the shock at radius Rq with Lorentz factor 7, and 
moving initially perpendicular to the shock, in the observer's frame; the shock has Lorentz 
factor Fs = ^s{R) ^ 1, so that most particles will be emitted perpendicular to the shock 
(Bednarz and Ostrowski 1998). We also have allowed for the shock deceleration. Because 
of a component of the magnetic field Bj_{R) = BiRi/R perpendicular to the shock normal, 
and thus to the particle initial direction of motion, the particle speed perpendicular to the 
shock Vz changes according to 

- « 1 - ^ rB^dR = 1 -^^ln{R/Ro) . (3) 

The shock moves instead with speed v shock /c ^ 1 — l/2r^. It is necessary that the shock 
overcomes the particle promptly; only in this case, in fact, the returning particle will see a 
shock with the same Lorentz factor as when it left the post-shock region, and conventional 
Fermi-Type acceleration will take place; if this were not the case, we would have a strongly 
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evolutionary shock, and the ordinary predictions of steady-state theory would not follow 
(and, most likely, the shock would not be an efficient particle accelerator). We must thus 
have VshockdR > v^dR, or, using the equations above, 

where R is the radius at which the particle recrosses the shock. Demanding again that the 
shock promptly overcomes the particle means that ln(i?/i?o) ~ 5R/Rq <ti 1, and thus we find 
that only particles with 

Ep = mTc" < 2TleB^R^ = E^^ (5) 

are accelerated by the usual Fermi mechanism. Particles with E > E^b will find a weaker, 
slower shock if they come back, and will not be significantly accelerated any more. Because 
of Eq. 10, we also see from the above that the condition that particles never come back 
(assuming that the PWB merges into the ISM field for it! > 10^^ cm) is only marginally 
stronger than the above: we find 

Ep ^ Emb = 2TleBiRi . (6) 

Gallant and Achterberg (1999) used R\ = Rd = 10^^ cm, the approximate radius at 
which the afterglow begins, the typical ISM field B — 1 /iG, and the initial shock Lorentz 
factor to obtain Eq. 2. But for a PWB, following the discussion above we may take for the 
field outside the fight cylinder radius (Eq. 31 of Goldreich and Julian 1969): 

.(«0..^io>'a.„^(i^)^in«M. p) 

Here P is the pulsar's rotation period, which we have scaled to 1 msec, which is about right 
for both neutron star mergers and SupraNovae; we took for the dipolar field of the neutron 
star a standard value, 6 is the beams' semiopening angle and '^{9) /I is a normalized fudge 
factor which Gooldreich and Julian (1969) do not estimate. However, inspection of Fig. 3 of 
Contopoulos, Kazanas and Fendt (1999) shows that '^{0)/I ^ (sin^)^; using 9 ^ 0.1 from 
Frail et al. (2001) and Panaitescu and Kumar (2002), we find 

B{Ri)Ri « 10^^ G cm , (8) 

independent of Ri. Strictly speaking, we have 100 > > 1, but for a Newtonian shock 
particle acceleration proceeds on a slower, diffusive timescale (Meli and Quenby 2002; see 
also the discussion later in this section). From Fig. 2 of Kirk et al. (2000), we see that the 
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particle spectral index departs from it asymptotic value for very large shock Lorentz factors 
when Fs 2. So we shall use 100 > > 1, but keeping this caveat in mind: 

10^^ eV<Emb<Qx 10^° eV , (9) 
depending on the instantaneous value of F. Since we know that (Piran 2000) 

r^ = r^(f)' (10) 

where F^ 100 (Panaitescu and Kumar 2002) is the initial shock Lorentz factor, we see that 
Ejnh spans continuously the range above, as R, the instantaneous shock position, changes. 

We thus have the following picture. So long as Emb > E^ax, particles cannot escape in 
the upstream direction: they are promptly overtaken by the shock because they are effectively 
kicked backwards by the environment magnetic field, and continue their acceleration cycle 
until they reach an energy E^ax- When this occurs, the particles cannot be effectively 
deflected by the ejecta's magnetic field: they will cross unhindered the hyperrelativistic 
shell, and be lost through the back of the shell. When this occurs, their energies in the 
laboratory (upstream) frame are a full factor 1/F^ lower than given by Eq. 1, which holds 
only for particles leaving the shell in the forward direction. However, when E^i, < Ej^ax-, 
particles which have been accelerated up to E^i, (in the laboratory=upstream frame) cannot 
be deflected backwards toward the shock any longer, and will be able to escape in the 
upstream direction without further acceleration. The critical moment occurs when E^nb — 
Emax (where now in Eq. 1 the shell Lorentz factor must be taken to vary as in Eq. 10): 
equating the two, we find that F 30, and thus that E^ax at this moment is a factor 1.5 
lower than Eq. 1 (and thus, within the accuracy of these estimates, it is identical to Eq. 1). 
We thus have that the energy range for particles escaping in the upstream direction is 

10^^ eV<E< E^ax ■ (11) 

It should be pointed out that the above picture differs from Gallant and Achterberg's (1999). 

We now have to derive the spectrum of these particles. Here there are two distinct 
possibihties. On the one hand, one may take Vietri's (1995) and Gallant and Achterberg's 
(1999) view that the particles which end up as UHECRs are the relativistic particles pre- 
existing the burst which must necessarily exist in a PWB. Under this assumption. Gallant 
and Achterberg showed that the spectrum will be of the form 

(h{E) oc E-^dE . (12) 

On the other hand, one may assume that the total number of energetic particles is 
dwarfed by the (formerly thermal) particles injected at the shock, in which case the spectrum 
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of particles can be derived as follows. Assume that a fraction q of the total fireball energy 
Ef is converted on the timescale tacc in particles of energy given by Eq. 6. There will then 
be a total amount of energy 

f^<t=^^ (13) 

in particles of energy ^ given by Eq. 6. Eliminating Tg and R, and using (Gallant and 
Achterberg 1999) 

(14) 



1 sGBc 

where is the upstream residence time, we find a particle spectrum 

dn (X E-^/'^dE , (15) 
Putting together all of the above, we shall take a spectrum of the form 

dnocE-^ ; 2 < 7 < 2.5 . (16) 

We shall then let 7 be a free a parameter in the range above, and use it to fit observations, 
in the next section. 

It is worth remarking at this point that the above results depends upon our having 
assumed that the GRB goes off inside a PWB (Konigl and Granot 2002), a natural assump- 
tion when GRBs are due to netron star mergers or to SupraNovae. While currently neutron 
star binary mergers are not favored as an explanation for long-duration GRBs, they appear 
attractive as a model for short GRBs; while SupraNovae are one of the leading contenders for 
the long bursts (see Dermcr 2002b for a review). But the above discussion applies, at least 
qualitatively as a lower limit, to a wide class of models, including Usov's (1992) model of a 
fast pulsar with exceptionaly strong magnetic field, models based upon ADIOS fiows (Bland- 
ford and Begelamn 1998,Blandford 2002), Spruit's models based on MHD winds (Drenkhahn 
and Spruit 2002, and references therein), or several distinct models all involving magnetized 
neutron stars (Duncan and Thompson 1992, Dai and Lu 1998, Blackman and Yi 1998). 

One may wonder whether losses in the somewhat more photon- and magnetic field- rich 
environment provided by a PWB make the estimate of Eq. 1 unrealistically optimistic. This 
was studied by Inoue, Guetta and Pacini (2002), who proved this to be the case only under 
exceptional circumstances. It was in fact their aim to explain the high-energy (> 1 GeV) 
emission occasionally seen by EGRET in GRBs, and this could be achieved only with rather 
extreme assumptions on the environment parameters. This is not at all in contradiction with 
the above discussion: GeV emission in GRBs is a rare phenomenon, so it makes sense that 
models to explain it require extreme values for the parameters in play. 
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A different critique was levied (Ostrowski and Bednarz 2002) concerning tfie spectral 
slope for particles accelerated at highly relativistic shocks. These authors pointed out that 
the claims of universality for the spectral index s ^ 4.2 (Kirk and Schneider 1987, Heavens 
and Drury 1988, Bednarz and Ostrowski 1998, Kirk et al., 2000) possibly overstate their 
case, because they neglect the fact that, as the regular magnetic field in the downstream 
region becomes stronger than the field associated with the turbulent component, the spectral 
slope must increase, the spectrum softens, and the claim to universality evaporates. This 
criticism, however, can be countered in two (orthogonal) ways. On the one hand, it is 
possible that it applies to a somewhat idealized situation, whereby post-shock material is 
incapable of generating a small-scale magnetic field, and thus the field compressed at the 
shock dominates the turbulent component. This is the view of Medvedev and Loeb (1999) or 
Thompson and Madau (2000), who proposed ways in which magnetic field of non-negligible 
amplitude could be generated around the shock, and would be intrinsically turbulent. On 
the other hand, one may take the view of Konigl and Granot (2002), according to whom 
the very large values of derived for some bursts (in particular, GRB 970508 has = 0.1 
as determined by several independent studies: Wijers and Galama 1999, Granot, Piran and 
Sari 1999, Ghevalier and Li 2000, Panaitescu and Kumar 2002) cannot possibly be accounted 
for in this way, but are instead easily accounted for by the compression of a magnetic field 
surrounding the burst explosion site, provided this occurs in a PWB, which again occurs 
naturally if bursts are due to binary neutron star mergers, or SupraNovae. The dominant 
component of this magnetic field is of course toroidal (Goldreich and Julian 1969), and its 
sweep up and compression would give a post-shock magnetic field mostly oriented along the 
shock surface, which would increase, not decrease, the shock effectiveness in the acceleration 
of particles. 

Another criticism was that the acceleration time scale was shorter than hypothesized 
in Vietri (1995), because the average energy gain < G > docs not scale as F^, where F 
100, Panaitescu and Kumar 2002) is the shock Lorentz factor. It was pointed out that 
this energy gain applies to the first shock crossing, but the following crossings would provide 
a more mundane < G 2 (Gallant and Achterberg 1999). This criticism is surely correct, 
but does not spoil the gist of the original argument, that the acceleration time scale is much 
shorter than for Newtonian shocks. Indeed, the fact that the energy gain per shock crossing 
is no longer infinitesimal, as it is for Newtonian shocks (where G — l = Vshock/c <^ 1), implies 
that the required number of shock crossings for acceleration to the EeV range becomes finite: 
from niT^c^ to the EeV range only ^ 20 shock crossings are required, in agreement with the 
cosiderable reduction of the acceleration time scale seen in numerical simulations (Meli and 
Quenby 2002). This very fact has led some of the severest critics of the scenario to propose 
themselves that UHEGRs are accelerated at highly relativistic shocks (Kirk, Achterberg, 
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Guthmann and Gallant 2000). 

There is a deep reason for this: we showed elsewhere (Vietri 2002) that, for particle 
acceleration around shocks of any speed, the following relation holds: 

P < G^-' >= 1 . (17) 

Here P is the probabihty of returning to the shock, G the energy amplification for one 
cycle across the shock, s the distribution function spectral index, and the average is to 
be taken after the raising to the 3 — s power; the above generalizes to relativistic shocks 
Bell's (1978) relation for the spectral index. The above relationship holds irrespective of the 
assumed scattering law. In the case of Newtonian shocks, both P and G differ from unity 
by quantities 0{v/c) ^ 1, and the two combine to give s = 4. In the hyperrelativistic case, 
both P and G differ significantly from unity, so that, since P ^ 1/2 (the post-shock fluid 
runs away from the shock with speed ~ c/3, so that it is easy for particles to be advected 
away), also G — 1 must grow away from its infinitesimal Newtonian values, and tends to 
G ^ 2. This, of course, shortens the acceleration time scale, and is a robust prediction of the 
theory of acceleration around hyperrelativistic shocks, independent of the scattering law. 



3. Energetics 

Another criticsm which has been levied against this scenario is that, contrary to our 
original claims, there is a mismatch between the observations of the energy injection rate 
in 7-ray photons due to GRBs, and in UHECRs. This claim has been made independently 
by Stecker (2000) and Berezinsky et al. (2002), who pointed out that the realization that 
the redshift distribution of GRBs should follow the star formation history in the Universe 
must lower the local energy injection rate. In particular, Stecker found the evolving ORB 
population to produce a 7 ray photon fiux which is an order of magnitude lower than that of 
just the fraction of UHECRs beyond 10^'^ eV, let alone the rest of the spectrum.. Waxman 
(2002) made some quite reasonable criticisms of the paper by Berezinsky et al. (2002). We 
consider here Stecker's work. We shall first show that Stecker's claims are not supported 
by currently available data, then derive the spectrum of UHECRs if these are emitted by 
GRBs, and then comment on the (rather large) uncertainties still surrounding the overall 
energetics of GRBs. 

Stecker's (2000) estimates are based upon old data: he used the AGASA fiux at > 
10^° eV (Takeda et al, 1998, Hayashida et al, 1999) which he took to be an unbroken 
power-law with index 2.78 beyond 10^° eV, while Abu-Zayyad et al (2002) explicitly state 
that their data beyond 10^^-^ eV depart from such a law, and state that the number of events 
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they observe is a full factor of 4/19.1 0.21 below the extrapolation of such law. Thus, 
while he took the local energy release rate in UHECRs to be 

HAGASA) = 2 X 10^^ erg Mpc'^ yr'^ , (18) 

HiRes data suggest a value a factor of 0.21 lower: 

eiHiRes) = 4 X 10^=^ erg Mpc'^ yr'' . (19) 

which compares already much better with the value he took for the energy release rate in 
7-ray photons, 2.2 x 10^^ erg Mpc^"^ yr~^. But even this value is superseded, since it is 
based on Schmidt's (1999) estimate of the GRB rate. In his new work, Schmidt (2001) gives 
the burst's rate (for isotropic bursts) as 

fiGRB = 0.5 Gpc-^ yr-' . (20) 

As correctly pointed out by Stecker (private communication), the above value includes short 
bursts, about which, strictly speaking, we know very little. In particular, we ignore what their 
overall energy budget is, and, most important to this paper, whether they have relativistic 
afterglows at all. So to play it safe, we reduce the above rate to exclude short bursts: 
According to Fishman and Meegan (1995) short bursts account for 30% of all bursts, and 
we consequently take the local rate of long GRBs as 

fiiGRB = 0.35 Gpc-^ yr-' (21) 

Schmidt does not give the average (isotropic) burst energy release, but we can get this from 
Frail et al. (2001), who, among many other things, give the isotropic burst energy release 
for all bursts with known redshifts. The reason why we prefer this over using Schmidt's 
work is that their total energy releases (taken from Bloom, Frail and Sari 2001) cover the 
spectral range 0.2 — 2QQQkeV , while Schmidt restricts himself to the range 50 — 300A;el^: 
there seems no obvious reason why we should adopt Schmidt's values, when more suitable 
ones are available. From Frail et al.'s Table I, we find 

Eav = 3.3 X 10^^ erg . (22) 

We thus find the energy release rate in 7-ray photons by GRBs as 

ecRB = 1.1 X 10^^ erg Mpc-^ yr-' . (23) 

Strictly similar results are obtained in an analogous way by using data from Panaitescu and 
Kumar (2000) and from Piran et al. (2001). Please notice that, while not all data from 
these three sets of authors are different (many bursts are in common) they fitted different 
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quantities with different methods, so the results are, for the most part, independent. This 

is larger than e^HiRes) by a hefty factor of ecRB I ^(mRes) = 3, and allows a full explanation 
of UHECRs by GRBs, even in the case in which the efficiency for particle acceleration is 
smaller than for photon production. 

We thus consider Stecker's (2000) objections superseded by new data. However, we 
are still left with the question of how much energy is produced by GRBs in the form of 
UHECRs in the whole energy range, Eq. 11, and not just for E > 10^*^ eV, and whether 
the observed spectrum of UHECRs at Earth can be fitted by a cosmological distribution 
of GRBs. In order to answer both of these questions, we repeat here a computation of the 
spectra of UHECRs as observed at Earth, for the evolving population discussed above. This 
does not mimick Waxman's (1995b) computations for several reasons. First, he uses a 
constant comoving density for GRBs, while, following Schmidt (2001), we use the second 
star formation history (SFR2) proposed by Porciani and Madau (2001); second, he uses a 
continuous energy loss approximation, while we take into account the discreteness of the 
photopion production events (but not of pair production); and third, what wc present here 
are simulations of observations (De Marco, Blasi, Olinto 2003) in order to stress that the 
fits to observations contain some amount of noise due to the finite number of events now 
available. For a full discussion of the method, simulations, and approximations, we refer 
the reader to De Marco, Blasi, Olinto (2003). Following Eq. 16,11, the particle spectrum is 
taken as 

n{E)dE ^E'-'dE ; E^^^ <E< E^^x (24) 

where 7 is a free parameter that is to be varied to allow fitting of Earth-based observations; 
Emin is from Eq. 11, while E^ax is irrelevant (provided it exceeds 10^° eV, which it does, 
Eqs. 1,11) because most energy is in the lowest energy particles (for 7 > 2). The spectrum 
we obtain is shown in Figs. 1 and 2 is obtained for the two values of 7 = 2.2 and 2.5. A 

smaller value of 7 yields only marginally different results. The steeper spectrum, like 7 = 2.5 
fits data below 10^^ eV much better, but, because of the obvious possibility of a Galactic 
contamination, we do not regard this accurate fit as either desirable or likely. We shall thus 
take the value 

7 = 2.2 (25) 

as our canonical result. For both values of 7 we under-reproduce the AG ASA counts beyond 
10^0 eV, but not the HiRes data. 

The local energy release rate in UHECRs, for the whole range of Eq. 11 is 

/ E \ 

euH = 1.5 X 10^5 — erg Mpc-' yr-' (26) 
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which exceeds Eq. 23 by a factor of 13. This overestimate is reduced when we consider 
the suggestion (De Marco, Blasi and Ohnto 2003) that the energies of AGASA events are 
overestimated by a full 15%, which brings their data in much better agreement with HiRes'. 
In this case, the total required energy release becomes 



which now exceeds Eq. 23 by a factor of 9 only. 

The estimate above is much lower than the claims by Scully and Stcckcr (2002), who 
hypothesized an excess by a factor of 10^~^, over the much smaller range E > 10^^ eV, while 
the above equation applies to E > 10^'' eV. The origin for the disagreement between our 
paper and Scully and Stecker's is twofold. First, we fit a different, flatter injection spectrum 
(2.2 rather than 2.75), but we still obtain excellent fits. Second, we use a €grb (Eq. 23) 
which differs from theirs. However, they do not compute it anew, but they refer to Stecker 
(2000). We have already shown that this estimate is based upon old data, and will thus not 
discuss their results any further. 

The slight overabundance of energy required to explain UHECRs as observed at Earth 
with respect to the energy injection rate in photons by GRBs, which we estimated above as 
included in the range 10 — 13, is not an embarrassment for the present model. Numerical 
simulations of radiative efficiencies at internal shocks, the currently favoured model for the 
prompt phase of GRB emission (Spada, Panaitescu and Meszaros 2000) show that these are 
~ 0.005, so that even the overproduction of UHECRs by a factor of 9 amounts to a total 
efficiency for UHECRs of a few percent, surely in line with conventional thinking (Draine and 
McKee 1993) about emission efficiencies for non-thermal particles around non-relativistic 
shocks. For relativistic shocks, given the shortening of the acceleration time-scale (Meli and 
Quenby 2002), acceleration should be, if anything, more efficient. 

For this reason, and for the reasons to be discussed in the next paragraph, we are not 
worried by this slight discrepancy. 



There are two important caveats which need to be made when discussing bursts' ener- 
getics. 




erg Mpc ^ yr 



-1 



(27) 



4. 



Two caveats 



Several authors have stressed that GRBs could have a redshift distribution strongly 
skewed toward large redshifts, but while very likely, this idea has not been confirmed yet. In 
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order to stress this point, we have taken the set of all GRB redshifts known to date, to see 
whether the low redshift distribution agrees better with a Porciani and Madau-like evolution 
law (Porciani and Madau 2001), or with a non-evolving one. We restrict our attention to 
small redshifts {z < 1) for two reasons: first, the Porciani and Madau evolution law is very 
distinctive at such low redshifts, and is not disputed by different authors; second, because 
it seems reasonable to assume that this sample is unbiased by the selection effects which 
obviously mar the higher end of the distribution. 

In Fig. 3 we show the observed cumulative distribution in redshift of known GRBs; data 
are taken from the compilation publicly available at http : / /www . mpe . mpg . de/ j cg/grbgen . html, 
which also gives the relevant references. We compare it with the theoretical expectation that 
the GRBs explosion rate follows the star formation rate. We consider the three different star 
formation histories SFl, SF2 and SF3 given by Porciani and Madau (2001). We use in the 
figure Qm — 0.3 and Qa — 0.7 and a Hubble constant Hq — BbhQ^km Mpc~^. For com- 
pleteness we compare the observed distribution also with a distribution cx (1 + 2;)^. A KS test 
shows that the observed distribution without GRB 980425 differs from the theoretical ones 
SFl, SF2, SF3 and oc (1 + z f with a significance of 24%, 25%, 28% and 32% respectively, 
while inclusion of GRB 980425 changes these values to 10%, 11%, 13% and 16%. 

Fig. 3 clearly shows that we have right now no compelling reason to believe that GRBs 
are characterized by very strong evolution with redshift, a uniform distribution providing an 
equally good, if not better, fit to current data. But furthermore, even if larger samples like 
the one to be collected by SWIFT should in the future provide a confutation of this (which 
we very much believe), still the point would remain that a perfectly normal subsample of 
a dozen or so objects may depart greatly from the average distribution, and appear like a 
uniform one. In this sense, then, the argument made by Stecker (2000) and Berezinsky et 
al. (2002) appears as an oversimplification of reality: in discussing sources of UHECRs, one 
cannot forget that the total number of sources one is dealing with is statistically very small, 
and likely to stay that way, and thus a proper discussion of the model on an energy basis 
should include the effects of fluctuations. 

Second, there is an important uncertainty in the global energetics of bursts, as we now 
discuss. Comparing the late-time evolution of bursts (which fixes the total energy in the 
kinetic form) with the energy released in 7-ray photons, it is possible to deduce the radiative 
efficiency for several individual bursts (Panaitescu and Kumar 2000). The radiative efficiency 
thusly derived is very high, always eceeding 66% (except for one case) and occasionally even 
approaching an incredible 90% (see Table III of Panaitescu and Kumar). We remark that 
these derived efficiencies always exceed the largest efficiency achievable in converting rest 
mass into photons around maximally rotating black holes (43%). Although this is not a 
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strict theoretical upper limit (models approaching 100% have been proposed, Lazzati et al, 
2000), still values derived from observations arc so large to throw doubts on the validity of 
the present version of internal shock model. In fact, the maximum efficiency that can be 
reached with this model is of the order of 20% (Guetta, Spada and Waxman 2001a) but only 
under ad hoc assumptions: for random distributions in the wind properties, the efficiency 
cannot exceed ~ .005 (Spada, Panaitescu, Meszaros 2000). Perfectly efficient models without 
internal shocks can be conceived (Lazzati, GhiseUini, Celotti and Rees 2000), but they cannot 
account for bursts' spectra (GhiseUini, Lazzati, Celotti, Rees 2000). 

For this reason, Kumar and Piran (2000) have proposed an ingenious modification of 
the fireball spectra, whereby hot spots in the emitting surface occasionally dominate the 
emission from the whole GRB. These hot spots are not typical of the whole surface; when 
our line of sight chances through one hot spot, we claim to see a burst, but when otherwise, 
we wouldn't recognize the event as such because too dim. In this way, they proposed that 
there be a significant number of objects (dissimulated GRBs) escaping our detection, and 
redressing the statistics. The amusing consequence of this is that, while of course this implies 
we would have to decrease the value of Eq. 23 (thus apparently going in Stecker's direction) 
there would be a substantially higher number of afterglows acting as acceleration sites. From 
Table II of Panaitescu and Kumar, we see that the apparently isotropic kinetic energy of 
bursts' afterglows is 

Eag = 2.5 X 10^^ erg . (28) 

This of course is strictly relevant to our acceleration mechanism, because as stated above, 
we are proposing that UHECRs are accelerated at external (= afterglow) shocks, so that 
this is the relevant source of energy to tap. Comparison with Eq. 22 shows indeed that the 
radiative efficiency, defined as Eav/ {Eav + Eag) ~ 0.6, very large indeed. Combining this 
with Schmidt's rate estimate, we see that there is a kinetic energy injection rate of 

iag = 1.3 X 10^^ erg Mpc^^ yr~^ , (29) 

comparable with Eq. 23. If the argument by Kumar and Piran (2000) is right, Eq. 29 is to 
be increased. In fact, the kinetic energy per burst in the afterglow is derived from very late 
observations, when the hot spots are not working any longer (Kumar and Piran 2000). hut 
the rate is to be increased to account for all bursts for which our line of sight docs not meet 
a hot spot. This factor is estimated by Kumar and Piran as roughly a factor of 100, bringing 
Eq. 29 in agreement with Eq. 27, with a vengeance. This incompleteness factor can be 
estimated also as follows: conventional models (Spada, Panaitescu and Meszaros 2000) find 
a radiative efficiency of ^ 0.005; observations give 0.6, so the hot spots bias our statistics 
by a factor of 100. 
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5. Doublets and triplet 

One argument often used against GRBs as sources of UHECRs (and, in general, against 
all explosive sources) is the presence of doublets and triplet in the data from AG ASA (Takeda 
et aL, 1999, Nagano and Watson 2000, Hayashida et al, 1997), where, occasionally, the lower 
energy event precedes the higher energy one. The argument goes that, since the time delay 
with respect to photons is a monotonically decreasing function of particle energy, no impul- 
sive source can have produced these clustered events, while obvioulsy no such problem exists 
with continuous sources. In the light of the present disagreement between the observational 
results of HiRes and AGASA, one may wonder whether the energy estimates even of the 
events in the clusters are correct. We shall not take this view, but accept the existence of 
the doublets and triplet as established, despite the caveat aJoove, and show that the purported 
inconsistency of the explosive models is not water tight. The argument in fact is an oversim- 
phfication, and neglects the fact that the time delay is a statistically distributed quantity, 
with both a mean and a variance, which we now compute to allow a further discussion. 

We take as a model for particle motion a Langevin-like equation 

dp dv V ^ 

-f. = m7— = e- A S (30 

dt dt c ^ ' 

where, that is, we assume 5 to be a Gaussian field with zero mean and known self-correlation 
function. Initially, we behave as if B were known exactly; we also use the assumption that 
the particle Larmor radius -Rl is large with respect to the typical distances over which the 
field is correlated, say Vc- Tc/Rl ^ 1, and return later to discuss this approximation. We 
take z as the axis of the initial motion, and decompose in approximate form the equations 
of motion as 

dvz Vy dv± 

777,7— ~ e—By - eB^— ; 7777-7- ~ ^B^V - ^ByX . (31) 
dt c c dt 

Here we have assumed that the particle's deflections from the unperturbed trajectory are 
small (which means that c,Vx ^ Vy <^ c) and thus the equations are valid to small- 

est significant order only; obviously, the fields are to be computed along the unperturbed 
trajectory, z. For a source located at distance D from the observer, we find 

Ay = I dZ I BJz)dz ; Ax = — I dZ I BJz)dz (32) 

^7cJo Jo rn-icjQ Jo 

and inserting this into the equation for the z component we find 

g2 rD i-z i-z 
/ dZ dz dy{B,{z)B,{z + y) + By{z)By{z + y))^ 

e2 rD 



/ {D-z)dz dy{B,{z)B,{z + y) + By{z)By{z + y)) (33) 
Jq Jo 
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where the last identity has been obtained by means of an integration by parts. We see that 
A|/, Ax, are hnear in the fields, and thus are distributed with zero mean, while the time 
delay (which is, by definition, the departure from an unperturbed motion with 7 3> 1, 
and thus equals the time delay with respect to photons), is quadratic in the fields and in 
the particle charge. Because of this, the time delay has a non-zero mean, which we now 
compute by means of the above equation: we obtain 

< Az >= --^ - z)dz r 2BUy)dy . (34) 

III' 1 C Jq Jq 

Here Bxx{s) —< Bx{x)Bx{x + sz) > is the rms (since we assumed it to have zero mean) of 
one component of the field pependicular to the line of separation, averaged over all field con- 
figurations; we also used isotropy to set B^^ = Byy, and homogeneity to drop the dependence 
of Bxx on X. We show in the Appendix that, under very reasonable physical assumptions, 

S,,(s) = ^So'5(-,«) (35) 

where 

^ 3rr»/fa)4/(™ + gf-igl) 

2/~!,-«/(!,)4/ 

For q; = 5/3 we would have the well-known Kolmogorov-Okubo law. 

When we insert the above into Eq. 34, we see that we have to compute the integral 



dygiy/rc, «) ■ (37) 

In the above, z is the (r^ cosmological) distance travelled by the particle; thus we expect 
z ^ Tc, since, as remarked in the Appendix, z is of the order of the radius of the Univers 
today, while Vc is of the order of the radius of the Universe at the time of the generation of 
the magnetic field. We can then approximate 

Jo Jo ^ ^Jo y fivjdy 8 

where we used Eq. 51. q{a) is illustrated in Fig. 3. Plugging the above into Eq. 34 we find 

q(a)e^B^D^rc , , 
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We are interested not just in the mean values, but also in the rms deviations from the 
mean. For the directions perpendicular to the z axis, we easily find with a computation 
analogous to the above 

V<(A.,>.<(A.,>^ ^^-^ . (40) 

Two particles arriving on the same experimental apparatus will, in average, be separated by 
\pi of the above, and this can be expressed in units of r^'- 

S„^,fM^(£Y'\ (41) 



This adimensional quantity Sn is the number of average correlation lengths which separate 
two average particles reaching the same experiment. It essentially measures whether two 
particles of the same energy, have travelled through the same patches of magnetic field: this 
is so, whenever 5n <^ 1, or otherwise when 5n>l. 

In the case 5n <^ 1, we may reasonably assume that the particles of the same energy 
emitted simulataneously in an explosive event, will arrive at the experiment also simultane- 
ously, with the time-delay derived above. In the case 5n ~ 1 or greater, since these particles 
will have travelled through patches of uncorrelated magnetic field, we cannot assume any 

longer that the particles will arrive simultaneously; a more reasonable assumption is that 
their time delays be distributed according to the statistics implicit in Eq. 33. We will present 
the full statistics elsewhere, and will limit ourselves here to derive the rms deviation of the 
time delay around its mean. We proceed as follows: using Eq. 33 we find 

iAzf = -^rdz r dz'{D-z){D-z') fds r ds' 
7 (' Jo Jo Jo Jo 

{By{z)By{z + s) + B,{z)B,{z + s)){By{z')By{z' + s') + B,{z')B,{z' + s')) . (42) 
We now average the above over all field configurations; to do so, we use Wick's theorem since 

— * 

we assumed S to be a Gaussian field: 

< {By{z)By{z + s) + B^{z)B^{z + s)){By{z')By{z' + s') + B^{z')B^{z' + s')) >= 
2{2B,,{s)B,,{s') + B,,{\z' - z\)B,,{\s - s'\) + B,,{\z - z! - ^\)B,,{\z + s - ^'|) . (43) 

We recognize in the first term to the right of the equahty the term which leads to (< /\z >)^; 
since it is our aim to derive the rms in the time delay, we drop it from now on. The second 
term on the rhs of the above can be seen to vanish (or more precisely, to give exponentially 
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small terms when Bu oc exp{—D/rc)), and will thus be neglected. Putting together all of 
the above, we find 

A_ ^< (A.)^ > -(< Az >r - \^^/^,^s (44) 
and we now define a convenient quantity as 

rms ^ 



<Az> V2 ' 



(45) 



which is the result we were searching for. It is worthwhile remarking that this is independent 
of all quantities, so long as our approximations (Eq. 38) are satisfied. 

We now may begin our discussion by remarking that propagation in the IGM clearly 
satisifies the hypothesis under which these results were derived, i.e., Tc/Rl <S 1 and homo- 
geneity, isotropy, and gaussianity (HIG, from now on) of the field (see Grasso and Rubinstein 
2001 for a review), so that the computation is immediately apphcable. We find for Eq. 41: 

5 -1 ^ ( ^ V U6) 

^ 2 X 10-11 G 100 Mpc 7 \l00rj ^ ' 

showing that, for very reasonable values of the parameters involved, particles which reach 
Earth experiments cross regions with uncorrelated magnetic field. In this case, these par- 
ticles' arrival times are spread out around their mean value (Eq. 33) with an rms of order 
A; if doublets and triplets include particles with energies different by a factor p, their mean 
arrival times differ by a factor 1 — 1/p^, and this corresponds to a splitting of (1 — l/p^)/A 
times the rms arrival time, for the lower energy ones. For typical observed values of p ~ 2, 
we find that the separation of the means in arrival times is only 3\/2/4 ?s 1, implying a 
separation of la only. This shows that there must be a superposition between the higher 
energy ones (assumed here for sake of argument to be travelling without dispersion in arrival 
times) and the lower energy ones (assumed instead to be spread out). Also, it is worthwhile 
to remark that the flux suppression at the lower energy due to statistics (~ e"^, as will be 
shown in aforthcoming paper) is nearly exactly compensated for by larger number of low en- 
ergy particles: for a spectrum going hke ^ at Earth, the product of the two, e~'^p^ ^ 2, 
again in agreement with observations. 

However, besides crossing the IGM, any particle reaching the Earth must also cross 
the Local Supercluster, and its associated magnetic field. Since the Local SuperCluster is 
obviously not a virialized structure, we may reasonably assume that HIG holds again. The 
current value of B is rather uncertain. If we take D = 10 Mpc, B 3 x 10~^ G and 
Tc — 1 Mpc as discussed by Blasi and Olinto (1999), we find rj Rl — 1/3, so that the above 
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computations are still (marginally) applicable; also 5n ~ 10. A more recent determination 
(Vallee 2002) has B ^ 1 fiG, and Tc = 100 kpc, so that again Vc/ Rl — 1, and the computation 
is marginally applicable; for these parameter values, 5n ~ 100. 

The bottom line of this argument is that, either in the IGM (for reasonable parameter 
values) or in the Local Supercluster (for currently favored, but still very uncertain param- 
eter values), the protons which reach the same experimental apparatus on Earth will have 
travelled through disconnected patches of magnetic field, and their time of arrivals will be 
spread around the mean by their large (Eq. 45) variance, and will mostly wash out the strict 
time ordering which is to be expected in the case of particles of (formally) infinite energy. 

A different, and important question, is whether the statistics implicit in Eq. 33 is 
consistent with the several instances of doublets, and triplet, observed, and the number of 
observed inversion in the ordering of arrivals. A full answer to this question requires working 
out the full statistics implicit in Eq. 33, and the comparison with the real data; the details 
of this work lie outside the scope of this work, and will be presented elsewhere. But it is 
nonetheless worth remarking that, on the basis of the arguments presented above, we expect 
roughly equal occurrences of multiplets with lower (higher) energy particles preceding higher 
(lower) energy particles. Prom Table 2 of Takeda et al. (1999) we see that this is precisely 
the case: in multiplets C1,C4,C5 the higher energy particle precedes the lower energy one, 
and the reverse occurs in multiplets C2 and C3. 



6. Comparison with other work 

Simultaneously with, and independently of Vietri (1995), two other papers (Milgrom 
and Usov 1995, Waxman 1995) proposed that UHECRs were accelerated in GRBs, but, 
while Milgrom and Usov did not propose a specific acceleration mechanism, Waxman (1995) 
proposed that Type II Fermi acceleration at the subrelativistic internal shocks could account 
for the acceleration of UHECRs. 

It is however still controversial whether the internal shocks can account for the GRB 
emission in the prompt phase. Spectra of GRBs display both a universal break feature at 
200 keV (Preece et al, 2000) and a wildly variable low energy spectral index (Crider et 
ai, 1997, Preece et ai, 1998, Ghirlanda, Celotti and Ghisellini 2002), which have proved 
difficult to reproduce theoretically: to wit, look at Fig. 1 of Meszaros and Rees (2000). In 
particular, it seems that pure synchrotron emission is incapable of explaining at least some 
(but not all! Tavani 1996, 1997) of the bursts' spectral properties (Ghisellini, Celotti, Lazzati 
2000). Possible alternatives have been discussed (Pilla and Loeb 1998, GhiseUini and Celotti, 
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1999, Panaitescu and Meszaros 2000) but none of them is widely credited with being able to 
explain the above-mentioned properties of the observed bursts. These difficulties should be 
contrasted with the success of the external shock model, which is widely credited with the 
explanation of the properties of the afterglow (Piran 2000). 

Since internal shocks are Newtonian (Waxman 2002), we do not expect the acceleration 
time scales to be short enough: the large average energy gain discussed above (G 2, Gallant 
and Achtcrbcrg 1999) only applies to relativistic shocks, and the ability of internal shocks 
to provide UHECRs in the short times allotted remains to be proved. This is especially 
so, since Waxman (1995, 2002) is appealing to the much less efficient Type II acceleration, 
rather than the standard Type I. But, even if Typell acceleration at a Newtonian shock 
should prove rather fast, still numerical simulations (see e.g. Meli and Quenby 2002) show 
explicitly that acceleration at relativistic shocks is much faster than Type I acceleration at 
Newtonian shocks, and this, in its own turn, is much faster than Type II acceleration (again 
at a Newtonian shock). Thus we may always expect Type I acceleration at external shocks 
to dominate Waxman 's assumed mechanism. 

There is no obvious reason to suspect that particle acceleration at internal shocks should 
be particularly effective: as emphasized by Waxman (2002), internal shocks are Newtonian, 
and thus a strong counterpressurc due to particles will result in a number of instabilities 
(Ryu, Kang and Jones 1993, and references therein) and in the smoothing out of the shock. 
The major consequence of this would be the suppression of the lower energy electrons which 
arc supposed to radiate away the winds' internal energy, and thus the muting of the burst. 
Acceleration of UHECRs at external shocks, instead, suffers no such problem ^: since it 
takes only a deflection by an angle C(l/r) to allow the shock to overtake a particle in the 
upstream region, as F — > oo the linear momentum communicated to the upstream section 
vanishes as l/F^, and no counterpressure/smoothing of the shock/instabilities ensue. 

The simple energy budget favors acceleration at external shocks. The apparently large 
radiative efficiency of bursts, contrary to all models for internal shocks (Spada, Panaitescu, 
Meszaros 2000), has prompted Kumar and Piran (2000) to propose an ingenious model to 
salvage internal shocks. According to them, hot spots in the shock surface strongly bias 
our detection abilities in their favor, and make it look like the overall radiative efficiency 
approaches unity; the many bursts without hot spots aiming at us, which are necessary to 
redress the statistic, would not be detected. In this model, the true radiative efficiency 
in the internal shock phase would be about 0.01 of what is observed, and there should be 
about 1/0.01 = 100 more GRBs than we detect. But what this means is that there would be 



We are indebted to P. Blasi for pointing this out to us. 
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100 times as many external shocks from which to draw UHECRs, because external shocks, 
for which beaming effects are smaller, (and thus also the hot spot mechanism would not 
be operational, Kumar and Piran 2000) suffer through no such criticism (Vietri 1998a). In 
other words: either one believes the largish 0.18 for GRB 970508) radiative efficiencies, 
in which case it is hard to believe in the internal shock model; or, one tries to salvage the 
internal shock model, in which case the energy balance is dominated by a large margin by 
the kinetic energy of the afterglow phase, which can be tapped by acceleration of UHECRs 
at external shocks. 

Apart from these objections, it is not easy to conceive at this point an observational 
test capable of distinguishing between the two models. The only hope appears to be the pro- 
duction of high energy neutrinos which must accompany the in situ acceleration of particles; 
occasionally, in fact, UHECRs will produce pions (and then neutrinos) through collisions 
with photons, in the moderately photon rich environment provided by the post-shock shells. 
If UHECRs are accelerated at internal shocks, the neutrinos thusly produced will arrive at 
Earth simultaneously with the photons of the burst proper and will have an energy ~ lO^^eV 
(Waxman and Bachall 1997, Guetta Spada Waxman 2001b). If UHECRs are accelerated at 
external shocks they will arrive at Earth simultaneously with the photons of the afterglow 
and will have a higher energy, > 10^^ (Vietri 1998a, 1998b). 



7. Conclusions 

In short, an analysis of the objections raised against acceleration of UHECRs at the 
external shocks of GRBs shows that none of them is, at the time of writing, dangerous for 
the model. 

The strongest objection levied so far is that due to Gallant and Achterberg (1999) that 
the highest energies of particles accelerated at external shocks is smaller than surmised in 
Vietri (1995); while in fact we pointed out that this objection is model-dependent, so is the 
reply, and ultimately the success of the scenario depends upon whether models in which a 
GRB explodes in a PWB (as discussed by Konigl and Granot 2002) or some magnetic-field- 
rich environment (Blandford and Begelman 1998, Usov 1992, Drenkhahn and Spruit 2002, 
Duncan and Thompson 1992, Dai and Lu 1998, Blackman and Yi 1998) are correct, or not. 

The strongest success of the model is once again the close agreement between the energy 
release in 7-ray photons and that necessary to account for UHECRs at Earth: to wit, 
compare Eqs. 27 and 29 (including the whole discussion, in Section 4). This agreement 
requires slightly different production efficiencies for photons and high energy particles, but 
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is otherwise unique to this model. 

It is also worth remarking that much information is contained in the observed properties 
of UHECRs: as noted above, in fact, multiplets contain information on the properties of the 
magnetic field along the line of sight which is peculiar to strongly time-varying sources, and 
which, coupled with other observed features such as the spectrum and the overall number of 
sources on the plane of the sky, can put the model through a severe test. A detailed analysis 
of this kind will be presented elsewhere. 

It is a pleasure to acknowledge extremely helpful conversations with P. Blasi. 



Appendix 

We wish to determine here the correlator for the Gaussian field: 

Bij{s) =< Bi{x)Bj{x + s) >■ . (47) 
We start out by considering the Fourier transform of Bij-. 

bij{k) = J B,^{s)e-'^-'dh . (48) 

Because of homogeneity, isotropy, and Maxwell's equation \j ■ B — Q one finds (Landau and 
Lifshitz 1987) 

bijik) = -^^{Sij - ^) (49) 

where E{k) is the mode energy distribution. If E{k) were to follow a pure Kolmogorov- 
Okubo law, we would have E{k) — Ak~^^^, but we have to modify this in two ways. First, 
the magnetic field outside of galaxy clusters cannot be due to the usual small scale turbulence 
that we are accustomed to in the laboratory, and thus the exponent in the above law need not 
equal 5/3. Second, if phase transitions at an early epoch (Grasso and Rubinstein 2001) are 
responsible for the generation of the magnetic field outside galaxy clusters, the field cannot 
be correlated over lengthscales exceeding the horizon scale at field generation. So we need 
to introduce a cutoff /(rg/c) for all modes with k < l/vc, for some Tc which we shall assume 
known from now on. We thus take 

E{k) = Ak-"f{rck) . (50) 

Strictly speaking, we do not know f{x) either; for computational purposes, in the following 
we shall take 

fix) ^ e-V- . (51) 
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Taking s — sz, we now compute B^xis) — Byy{s). Prom the inverse of Eq. 48 we find 

BUs) = j (52) 
which, using Eqs. 50 and 49 yields 

Bxx{s) = ^ / k-'^f{rck){cos^ 9 + sin^ ^sin^ 0)6^'=^"°^^* dcf) dcosO = (53) 

^ rk-^f{nk)dk f ^ + '°''^ '^--^dcosg. 

2 Jo J-i 2 



We find 

B 



/V 9 . /"°° , ^ , V , /sin /cs cos ks sin /cs\ . 

(.) = j . (54) 

When s — s> 0, because of isotropy obviously B^x -Sq/S, where Bq is the rms field. Letting 
s — > in the above equation, we find 

/■oo 

B^ = 27r^A / k-'^f{rck)dk (55) 
Jo 



which allows us to rewrite 



5..(«) = gi^o 9{-, a) ; ,(., a) = 2j-y-^fiy)dy ^''^ 



and similarly for B, 



yy 
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Fig. 1. — Comparison of expected and observed spectra of UHECRs at Earth. Solid line, 
analytical computation of the spectrum. Crosses with bars, numerical simulations of the 
AGASA data. Circles with bars, AGASA data. Suqares with bars, HiRes data. The bursts' 
redshift distribution is Porciani and Madau's (2001) SFR2. The injection spectrum extends 
over the range given by Eq. 11; the injection spectral slope is 7 = 2.2. 



'lo- 



ts 10" 



10" 



+ ^ t__ 



AGASA 


Y=2.5 


E>10" 


866 


E>10'"' 


55.38 ±7.57 


E>10'" 


2.54 ± 1.68 



18.5 



19 



19.5 



20 



20.5 21 
loglO(E) [eV] 



Fig. 2. — Same as Fig.l, except that 7 = 2.5. 
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Fig. 3. — The observed cumulative distribution function in redhsift, for all GRBs with z < 1. 
The two observed samples exclude (Sample 1) or include (Sample 2) the low redshift burst 
GRB 980425. These observed distributions are compared with theoretical distributions, 
obtained for three different star formation histories, taken from Porciani and Madau 2001: 
SFl (continuous curve), SF2 (long-dashed curve), SF3 (dashed curve). The dotted curve 
shows a non-evolving distribution oc (1 + 2;)^. 



